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1. Introduction 

In this brief review, I summarize progress in recent years in the thermodynamics 
and collective excitations of quantum chromodynamics (QCD) based on the hard- 
thermal-loop (HTL) effective theory. This is a very active research area, especially 
due to the pending challenges from the ultra-relativistic heavy-ion collision experi¬ 
ments at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Labo¬ 
ratory, the Large Hadron Collider (LHC) at the European Organization for Nuclear 
Research, as well as the forthcoming experiments and the Facility for Antiproton 
and Ion Research (FAIR) at the GSI Helmholtz Centre for Heavy Ion Research. 

Asymptotic freedom^’^ predicts that at sufficiently high temperature (T) and 
density/chemical potential (/i) hadronic matter would undergo a phase transition 
(or a crossover) to a novel deconhned phase of quark-gluon plasma (QGP). The 
launch of heavy-ion experiments at the RHIC and most recently at the LHC opened 
an unprecedented era for the study of matter under extreme conditions: the initial 
temperature at the RHIC is expected to reach up to 2Tc, while that of the LHC 
up to dTc, where ~ 160 MeV is the pseudo-critical temperature for the QCD de¬ 
confinement transition. Strikingly, the matter created in the collisions behaves like 
a nearly ideal liquid characterized by a small shear viscosity, which has spoiled the 
naive expectation from asymptotic freedom for a weakly coupled QGP (see Ref. 3 for 
a review). The QCD running coupling expected in the experimental energies is on 
the order of unity g ~ 0(1), which is some intermediate value neither infinitesimally 
small nor infinitely large. There have been insights gained from the strong-coupling 
formalism based on the anti-de Sitter/conformal field theory (AdS/CFT) corre¬ 
spondence (see Ref. 4 for a review), however in order to have a direct access to 
the physical mechanisms of the QGP thermodynamics and realtime dynamics and 
considering the fact that finite density and realtime dynamics are yet to be settled 
challenges in lattice QCD (see Refs. 5,6 for reviews), continuum QCD methods such 
as resummed perturbation theory are indispensable. 

Collective excitations in the QGP introduce daunting challenges in practice. In 
addition to the intrinsic energy scale T, collective behaviors of the QGP generate 
two thermal scales, namely the (chromo)electric scale gT and the (chromo)magnetic 
scale g^T. The electric scale is dominating at high T where g is small and it is well 
described by the hard-thermal-loop effective theory^^^^ (see also Refs. 13-19 for 
reviews on various aspects of the developments). The magnetic scale, which con¬ 
tributes more and more significantly as T is decreased towards the phenomenologi¬ 
cally relevant regime where g ^ 0{1), is beyond the scope of conventional resummed 
perturbation theory due to the so-called Linde problem. por heavy-ion experi¬ 
ments to have the greatest possible impact on science, it is essential to make as close 
as possible a connection to the fundamental theory QGD. There is thus an urgent 
need for theoretical frameworks based rigorously on QCD which can be applied to 
the QGP in the phenomenologically relevant regime. In the following, I briefly re¬ 
view the recent developments for this concern from both the electric and magnetic 
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scales based on the HTL effective theory. In Sec. 2, I discuss the progress in the 
electric scale by a next-to-next-to-leading order calculation of QCD thermodynam¬ 
ics at finite T and ^ from the hard-thermal-loop perturbation theory. In Sec. 3, I 
discuss the progress in the magnetic scale by a first HTL analysis of the collective 
excitations of hot QCD including the magnetic scale. I conclude in Sec. 4 with an 
outlook for future perspectives. 


2. Progress in Electric Scale: Thermodynamics 

It had been a longstanding problem that the resulting series of a weak-coupling 
expansion is poorly convergent at finite T unless the coupling is tiny,^^“^^ thus 
a straightforward perturbative expansion in powers of g for QCD is not of any 
quantitative use for the temperatures achieved in heavy-ion collisions. The poor 
convergence stems from the fact that at high T the classical solution is not well- 
described by massless degrees of freedom, and it is instead better described by 
massive quasiparticles with non-trivial dispersion relations and interactions gener¬ 
ated by the thermal scales. This calls the need for reorganizing the perturbative 
series which treats the thermal scales more carefully. Hard-thermal-loop perturba¬ 
tion theory (HTLpt) is a reorganization scheme for thermal QCD incorporating the 
electric scale.It systematically shifts the expansion to being around an ideal gas of 
quasiparticles with a thermal mass on the order of gT. HTLpt is a gauge-invariant 
generalization of screened perturbation theory which is a reorganization scheme 
for scalar field theory at finite T inspired in part by variational perturbation the¬ 
ory The thermodynamic calculations from HTLpt have been firstly carried out 
at finite T and vanishing p, at one-loop or leading order two-loop or 

next-to-leading order (NLO),^®’^^ and three-loop or next-to-next-to-leading order 
(NNLO).^®^^^ In recent years, the corresponding calculations have been generalized 
to finite T and The HTLpt framework has been applied to evaluating the 

equation of state for cold dense quark matter (i.e. vanishing T and finite with 
application to compact stars.Application of some HTL motivated approaches to 
thermodynamics and various susceptibilities can be found in Refs. 63-72. In the fol¬ 
lowing of this section, I briefly discuss the setup of HTLpt and the recently obtained 
NNLO thermodynamic potential at finite T and p from which various thermody¬ 
namic functions and susceptibilities are derived.®® 


2.1. Hard-Thermal-Loop Perturbation Theory 

Hard-thermal-loop perturbation theory is a reorganization of the perturbative series 
of thermal gauge theories. The HTLpt Lagrangian density for QCD in Minkowski 
space can be written as 


C = (/IqCD + ^HTL)|g^yjg + AThtl • 


(1) 
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Here >Cqcd is the QCD Lagrangian density that reads 

>Cqcd = — + >Cgh + >Cgf + A£qcd , (2) 

where G^'' = A'^ — d''A^^— ig[A>^, A'^] is the gluon field strength tensor, = d^ — 

igA^ is the covariant derivative, and the term with the quark fields ip contains an 
implicit sum over the Nf quark flavors. The term A/Iqcd contains the counterterms 
necessary to cancel ultraviolet divergences in perturbative calculations. The ghost 
term Cgh depends on the form of the gauge-fixing term Cg{. 

The gauge-invariant HTL improvement term /Ihtl reads^^ 

CuTh = (1 - (^%). , 

( 3 ) 

where = (l,y) is a light-like four-vector with y being a three-dimensional unit 
vector and the angular bracket indicates an average over the direction of y. The two 
parameters mo and niq can be identified with the Debye screening mass and the 
quark thermal mass, respectively, and account for screening effects. HTLpt is defined 
by treating 5 as a formal expansion parameter. By coupling the HTL improvement 
term (3) to the QCD Lagrangian (2), HTLpt systematically shifts the perturbative 
expansion from being around an ideal gas of massless particles to being around an 
ideal gas of massive quasiparticles which are more appropriate physical degrees of 
freedom at high temperature and/or density. 

Physical observables are calculated in HTLpt by expanding in powers of 5, trun¬ 
cating at some specified order, and then setting 5 = 1. This defines a reorganization 
of the perturbative series in which the effects of mjj and terms in (3) are in¬ 
cluded to leading order but then systematically subtracted out at higher orders in 
perturbation theory by the Smjy and terms in (3). The HTLpt Lagrangian (1) 
reduces to the QCD Lagrangian (2) if we set 5 = 1. If the expansion in 6 could 
be calculated to all orders, the final result would not depend on mo and rriq when 
we set 5=1. However, any truncation of the expansion in 5 produces results that 
depend on mo and iriq. As a consequence, a prescription is required to determine 
mo and vriq as a function of T, g and as- Note that HTLpt is gauge invariant 
order-by-order in the 5 expansion. 

The HTLpt expansion generates additional ultraviolet divergences. In QCD per¬ 
turbation theory, renormalizability constrains the ultraviolet divergences to have a 
form that can be cancelled by the counterterm Lagrangian A£qcd- There is yet 
no proof for the renormalizability of the HTL perturbation expansion, it has been 
shown in Refs. 48-52,59,60 that it is possible to renormalize the HTLpt thermody¬ 
namic potential through NNLO with a counterterm Lagrangian A/1htl containing 
only a vacuum counterterm, a Debye mass counterterm, a fermion mass countert¬ 
erm, and a coupling constant counterterm. The necessary counterterms for renor- 
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malization of the NNLO thermodynamic potential are 


Afo 

Am^ 

Arrig 

SAas 


q ; s ( 5(1 — 5)mjj 


<:\2 4 

llCA — 4 sf 
127re 

- 6)ml 

Sttc CA 

IIca ~ 4sf 2 r2 


(4) 

(5) 

( 6 ) 
(7) 


where ca = Nc, cIa = — 1, sf = A//2, dp = NcNf, and S 2 f = CpSf with 

Cp = {N^ — l)/2Nc. Note that the coupling constant counterterm (7) is consistent 
with one-loop running of as- 


2.2. NNLO Thermodynamic Potential at Finite T and /jl 

The calculation of the thermodynamic potential in HTLpt involves the evaluation 
of vacuum diagrams (see Figs. 2 and 3 in Ref. 52 for the diagrams through NNLO). 
The fact that mp and to, are of the order gT suggests that mp/T and mq/T can 
be treated as expansion parameters of order In order to make the calculation 
of the thermodynamic potential analytically tractable in practice, a mass expansion 
in terms of mp/T and mq/T is carried out after the 5 expansion is done. It was 
shown that the first few terms in the mass expansion gave a surprisingly accurate 
approximation to the exact result. 45 resulting NNLO thermodynamic 
potential is completely analytic, and it is accurate to order g^ in the weak-coupling 
limit. Defining H(z) = 4'(z) -|- 4'(z*) with z = 1/2 — ifi and 44 being the digamma 
function, Og = —dA'i^'^T^/rh^ and x = x/2'kT for dimensionless variables, the NNLO 
thermodynamic potential for QCD at finite T and y (with /ij being the chemical 
potential for quarks with flavor /) reads®® 
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where the sums over / and g include all quark flavors and H^nlo pure-glue 

contribution to the thermodynamic potential that reads^® 


^NNLO _ 1 I 
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(9) 


As discussed in Ref. 52, the two-loop perturbative electric mass for gluons in¬ 
troduced by Braaten and Nieto in Refs. 29,30 originally at finite T and vanishing /r 
is the most suitable one for NNLO HTLpt calculations, and it will be thus adopted 
in the next sections. The finite T and /x generalization was obtained in Ref. 73 and 
the resulting mjj reads 
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( 10 ) 


The effect of the in-medium quark mass parameter niq in thermodynamic functions 
is small and it is thus set to 0 which is the three-loop variational solution for vriq 
following Ref. 52. 


2.3. Thermodynamic Functions 

In this section I review results for the NNLO HTLpt pressure, trace anomaly, and 
speed of sound obtained in Refs. 52, 60. In all the results, the one-loop running 
corresponding to Eq. 7 of the coupling renormalization is used. It reads 

«.(A) = ^. (11) 

bot 

with t = ln(A^/A/^) and bo = (IIca — 2A^/)/127r. The scale A^jg is fixed by 
requiring that as(1.5GeV) = 0.326 which is obtained from lattice measurements.^^ 
This gives Aj^ = 176 MeV. For the renormalization scale, Ag and A^ are used 
for purely-gluonic and fermionic graphs, respectively. The central values of these 
renormalization scales are taken to be Ag = 27rT and A = Aq = 27r-y/T^ -k/x2/7rT 
In all plots the thick lines indicate the result obtained using these central values and 
the light-blue bands indicate the variation of the results under variation of both of 
these scales by a factor of 2, e.g. ttT < Ag < AttT. For all numerical results below 
Ca = Afc = 3 and Nf = 3 are set. 
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Fig. 1. Comparison of the Nf = 2 + 1, = 0 (left) and = 400 MeV (right) NNLO HTLpt 

pressure with lattice data from Borsanyi et and Bazavov et al.^^ 


2.3.1. Pressure 

The pressure is a key quantity from which all other thermodynamic functions 
can be derived. It is obtained directly from the thermodynamic potential by 

V = —nNNLO ■ (12) 

Fig. 1 shows the scaled NNLO HTLpt pressure for = 0 (left) and = 400 MeV 
(right) with lattice data from Refs. 75-77. As we can see from this figure, there is 
quite good agreement between the NNLO HTLpt pressures and the lattice data for 
T > 200 MeV when the central value of the scale is used. Since HTLpt does not 
incorporate the center symmetry Z(Nc), there is no reason to expect agreement 
with the lattice data at temperatures close to and the agreement seen may be 
fortuitous. 


2.3.2. Trace Anomaly 

The trace anomaly is obtained from the pressure by 

This quantity is the trace of the energy-momentum tensor and vanishes for an ideal 
gas of massless particles, it thus measures the breaking of conformal symmetry by 
quantum effects for such a system. Fig. 2 shows the NNLO HTLpt trace anomaly 
scaled by for = 0 (left) and fiB = 400 MeV (right) together with lattice data 
from Refs. 76, 77. As we can see from this figure, there is quite good agreement 
between the NNLO HTLpt trace anomalies and the lattice data for T > 220 MeV 
when the central value of the scale is used. 
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Fig. 2. Comparison of the Nf = 2 + 1, = 0 (left) and pB = 400 MeV (right) NNLO HTLpt 

trace anomaly with lattice data. The pQ = 0 lattice data are from Ref. 76 and the pB = 400 MeV 
lattice data are from Ref. 77. 




Fig. 3. Comparison of the Nf = 2 + 1, pB = 0 (left) and pB = 400 MeV (right) NNLO HTLpt 
speed of sound with lattice data. The pB = 0 lattice data are from Ref. 76 and the pB = 400 MeV 
lattice data are from Ref. 77. 


2.3.3. Speed of Sound 

The speed of sound Cg is another phenomenologically relevant quantity defined as 


c 


2 

s 


(fp_ 

~de' 


(14) 


Fig. 3 shows the NNLO HTLpt speed of sound for = 0 (left) and = 400 MeV 
(right) together with lattice data from Refs. 76,77. As we can see from this figure, 
there is quite good agreement between the NNLO HTLpt speeds of sound and the 
lattice data when the central value of the scale is used. 
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Fig. 4. Comparison of the Nf = 2 + 1 NNLO HTLpt results of the 4th order diagonal single 
quark number susceptibility (left) and the only non-vanishing 4th order off-diagonal quark number 
susceptibility (right) with lattice data. In the left figure the dashed blue line indicates the Stefan- 
Boltzmann limit for this quantity. The data labeled BNL-BI(uudd), BNL-BI(u,s), BNL-BI(uuss), 
and TIFR come from Refs. 78, 79, 80, and 81, respectively. 

2.4. Susceptibilities 

The quark number susceptibilities are another set of phenomenologically relevant 
quantities. These functions carry information about the response of the system to 
nonzero density. By taking derivatives of the pressure with respect to chemical 
potentials, we obtain the quark number susceptibilities 


where p, = {pu, Pd, ■■■, PNf) representing a separate chemical potential for each 
quark flavor. The left panel of Fig. 4 shows the 4th order single quark susceptibility 
^uuuu eomparing to lattice data from Refs. 78-81. As can be seen from this figure, 
the scale variation of the HTLpt result is quite small for this particular quantity and 
there is good agreement with the lattice data. The right panel of Fig. 4 shows the 
4th order off-diagonal single quark susceptibility which is also in reasonably 

good agreement with the lattice data. 

Since the directly accessible information in the experiments are baryon num¬ 
ber fluctuations rather than quark number fluctuations, we should pay particular 
attention to the baryon number susceptibilities defined as 



(16) 


with pb = Pu + Pd + Ps- The left panel of Fig. 5 shows the scaled 2nd order baryon 
number susceptibility comparing to lattice data from Refs. 78,79,81-83. As can been 
seen from this figure, the NNLO HTLpt result is in good agreement with the lattice 
data for T > 300 MeV. The right panel of Fig. 5 shows the scaled 4th order baryon 
number susceptibility comparing to the lattice data. The NNLO HTLpt result is 
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Fig. 5. The scaled 2nd (left) and 4th (right) order baryon number susceptibilities compared with 
various lattice data. The lattice data labeled WB, BNL-BI(B), BNL-BI(u,s), MILC, and TIFR 
come from Refs. 82, 78, 79, 83, and 81, respectively. 


consistent with the lattice data shown, however the lattice error bars on this quantity 
are somewhat large and the data are restricted to T < 400 MeV, making it difficult 
to draw firm conclusions from this comparison. That being said, HTLpt makes a 
clear prediction for the T-dependence of the 4th order baryon number susceptibility. 
It will be very interesting to see if future lattice data agree with this prediction. 

3. Progress in Magnetic Scale: Collective Excitations 

As mentioned in Sec. 1, conventional thermal perturbation theory breaks down 
at the magnetic scale g^T due to the Linde problem.The nonperturbative 
nature of the magnetic scale is intimately related to the confining property of the 
dimensionally reduced Yang-Mills theory at high temperature. This suggests that 
a confinement mechanism should be incorporated within perturbative expansions 
even when dealing with the deconfined QGP phase. There have been many efforts 
over the last two decades trying to reconcile resummed perturbation theory with 
the magnetic scale and this is still a key open question in the field.Color 
confinement is deeply related to positivity violation of the spectral function: if the 
spectral function of a particle is not positive semi-definite, no Kallen-Lehmann 
representation exists, it is then not part of the physical spectrum and thus confined 
(see Ref. 91 for a review). This problem has been studied intensively using lattice 
QCD and functional methods at both zero and finite temperatures for gluons (see 
Refs. 92-100 and references therein, see also Refs. 101,102 for reviews). The study 
of the quark sector has not been equally conclusive. 

Conventional thermal field approaches to hot QCD are based on massive quasi¬ 
particles which only generate short-range correlations. In order to describe a 
strongly coupled QCP, long-range correlations, whose carriers are light and/or mass¬ 
less modes, are a crucial ingredient. There have been hints on massless modes in a 
QCP from functional methods.At the thermal held frontier, there has been a 
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series of studies on massless modes in Nambu-Jona-Lasinio model, Yukawa model, 
QED and In the following, I briefly review the first study on massless 

modes in hot QCD using confining gluons recently reported in Ref. 114 which shows 
genuine non-Abelian features such as positivity violation. 


3.1. Gribov-Zwanziger Formalism at Finite T 


A formalism to tackle the issue is the Gribov-Zwanziger (GZ) action, which is well- 
known from the study of color confinement. regulates the IR behavior of 
QGD by fixing the residual gauge transformations, i.e., Gribov copies (see Ref. 117 
for a review), that remain after applying the Faddeev-Popov procedure. The GZ ac¬ 
tion is renormalizable, and it thus provides a systematic framework for perturbative 
calculations (i.e., (/ <C 1) incorporating confinement effects. The gluon propagator 
in general covariant gauge reads 





( 1-0 


pr pv 
P2 


p2 

P^ + lt ’ 


(17) 


where ^ is the gauge parameter and the Landau gauge ^ = 0 has been well explored 
in practice (see Refs. 118,119 for reviews). The Gribov parameter 7 g is solved self- 
consistently from a gap equation that is defined to infinite loop orders. The GZ gluon 
propagator is IR suppressed, manifesting confinement effects, and it is a signihcant 
improvement over the one from the Faddeev-Popov quantization which forms the 
basis for conventional perturbative calculations. The gap equation at one-loop order 
can be solved analytically at asymptotically high T and gives^^°’ 


7g = 


D-1 N, 2 
D 4v^7r^ 


(18) 


where D is the space-time dimensions. Eq. (18) provides a fundamental IR cutoff at 
the magnetic scale for the finite-T GZ action. The effectiveness of the GZ framework 
in the study of Yang-Mills thermodynamics has been explored in Refs. 120-123. 


3.2. Quark Thermal Self-Energy 

An important measure for the collective behavior of a QGP is the self-energy of 
quarks and gluons, from which thermal masses, dispersion relations, and spectral 
functions of collective excitations are derived. The Euclidean one-loop quark self¬ 
energy reads 

^{P) = {lyfCpi rS{K)YD^’'{P-K) , (19) 

1{K} 

where S{P) = l/f* is the quark propagator, and D^^''{P) is the gluon propagator 
which is taken from Eq. (17). It is worth noting that there have been similar studies 
for the quark self-energy with nonperturbative gluons at finite density^^^’ and in 
strong magnetic fields. 
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Fig. 6. The quark thermal mass mq^'yc) from Eq. (21) scaled by the perturbative value mg(0). 


When g < 1 , the leading contribution from E(P) to the two-point correlation 
function is from soft external momenta P ~ gT, and the leading contribution to the 
loop integral in Eq. (19) is from k ~ T.® This suggests that for studying the high-T 
behavior of the self-energy in the small-coupling regime, we may expand Eq. (19) in 
terms of small P following the systematics of the HTL effective theory. As a result, 
the gauge-invariant contribution to Eq. (19) reads^^^ 


S(P) 


_ pOO 

-{igfCp^ i 

do 



dfl h±(fc,7G) 

4E5, 


X 


iyo -I- fc • 7 


iPo + k — ^ 


-h 


i'jo - k-j 


iPo — k + ^ 


( 20 ) 


where k = k/k with k = |fc|, E^ = y/k"^ ± 17 ^, h±{k, 7 g) = n_B(\/fc^ ± ij'^)+nFik) 
with ns and np the Bose-Einstein and Fermi-Dirac distributions, and f dO = 
d(()/J'dcos 6 ». 


3.3. Quark Thermal Mass 


The quark thermal mass incorporating effects from the magnetic scale reads 


w^(7g) 


g^Cp 


E 

± • 


dk 


k'^h±{k,-fG) 




( 21 ) 


which reduces to the conventional HTL one, TOg(O) = Cpg‘^T^/8, for 7 g = 0. The 
scaled quark thermal mass raq(^Q)/mq{Q) is shown in Fig. 6 . It is clear from the 
figure that mq receives negative contributions from 7 g, which is a manifestation 
of anti-screening effects generated by the magnetic scale. Although the effect is 
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Fig. 7. Dispersion relations (upper panel) and the corresponding residues (lower panel) for the 
particle (uj+), plasmino (a;_) and Gribov (uig) poles. 


modest in the studied range of couplings, this is a profound signal of the build-up 
of long-range correlations in the system and similar anti-screening effects have been 
observed on the lattice for the Debye screening mass.^^^ 


3.4. Massless Mode and Positivity Violation 

The dispersion relation is obtained by analytically continuing the self-energy 
Eq. (20) to Minkowski space and then solving the poles in the corresponding quark 
propagator iS~^{P) = p—T,{P) = 0. The resulting dispersion relations and residues 
of the poles are displayed in the upper and lower panels of Eig. 7. 

In contrast to the conventional HTL case, there are three poles in the propagator. 
Eirstly, the screened quasi-particle excitations are recovered, 

‘^ = w+(p;7g), w = w_(p;7g), (22) 

the so-called particle and plasmino a;_ modes, with a;±( 0 ; 7 G) = niq(^a) as 
expected. Both uj±/mg{'yc) and Z± are ^-independent, and this has been verified 
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explicitly up to ~ 2 in Ref. 114. This property is exactly the same as in the 
conventional HTL effective theory, and it is thus a non-trivial consistency check of 
the setup. 

In addition to the massive modes, there exists a novel excitation named Gribov 
pole as in Ref. 114, 


w = ujg{p-,^g) ■ 


(23) 


It describes massless fermionic excitations in the plasma with dispersion relation 
u! = VsP at small momenta, with Vg ~ l/-\/3 (speed of sound) independent of 
g for the studied range. The Gribov mode “grows” in the (a;,p)-plane while the 
magnetic scale is increasing (through increasing g), and this effectively introduces 
a new magnetie scaling behavior to the non-Abelian plasma. The vertical lines in 
Fig. 7 schematically demonstrate how the Gribov mode grows: at small coupling, 
e.g., g = 0.5, the mode terminates at rather small momentum; as the coupling 
increases, to e.g., g = 2, the permitted momentum range increases accordingly. At 
larger momenta than the permitted ones for each coupling, we are hitting branch 
cuts and Landau damping takes place as a consequence. The Gribov pole goes 
along with a residue Zg{p) < 0 which directly implies positivity violation of the 
corresponding spectral functions in the region of space-like momenta. These novel 
features are direct manifestations of long-range confinement effects surviving at 
finite T in the non-Abelian plasma. The results reflect common features of Gribov- 
like approaches,though the calculation was done via the GZ action. 


4. Conclusions and Outlook 

In this brief review, I have attempted to discuss recent progress of hard-thermal- 
loop effective theory on both the electric and magnetic scales of a hot QGD plasma. 
The HTLpt thermodynamics calculation has been a daunting task, and the NNLO 
thermodynamic potential at finite T and pt resembles a continuous effort over the 
past 15 years. Much confidence has been gained from the results that HTLpt may 
provide a good description for QGP thermodynamic functions and various suscepti¬ 
bilities at T > 300 MeV, and it would be interesting to apply the HTLpt framework 
to realtime quantities at these temperatures. 

Comparing to the electric sector, the magnetic sector has been explored to much 
less degree. The uncovering of the massless Gribov mode has been an exciting 
attempt in exploring the significance of the magnetic scale to a non-Abelian plasma. 
It sheds new light on the active degrees of freedom released in course of a heavy- 
ion collision through which a strongly coupled QGP might emerge. It would be 
extremely tempting and challenging to explore the phenomenological significance of 
this new mode in interpreting experimental data. 
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